Bose-Einstein Condensation of Pions in High Multiplicity Events 
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We present microcanonical ensemble calculations of particle number fluctuations in the ideal pion 
gas approaching Bose-Einstein condensation. In the samples of events with a fixed number of all 
pions, N-k, one may observe a prominent signal. When N n increases the scaled variances for particle 
number fluctuations of both neutral and charged pions increase dramatically in the vicinity of the 
Bose-Einstein condensation line. As an example, the estimates are presented for p + p collisions at 
the beam energy of 70 GeV. 
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The phenomenon of Bose-Einstein condensation 
(BEC) was predicted long time ago Tremendous 
efforts were required to confirm BEC experimentally. 
The atomic gases are transformed into a liquid or solid 
before reaching the BEC point. The only way to avoid 
this is to consider extremely low densities. At these 
conditions the thermal equilibrium in the atomic gas 
is reached much faster than the chemical equilibrium. 
The life time of the metastable gas phase is stretched 
to seconds or minutes. This is enough to observe 
the BEC signatures. Small density leads, however, to 
small temperature of BEC. Only in 1995 two exper- 
imental groups succeeded to create the 'genuine' BE 
condensate by using new developments in cooling and 
trapping techniques 0. Leaders of these two groups, 
Cornell, Wieman, and Ketterle, won the 2001 Nobel 
Prize for this achievement. 

Pions are spin-zero mesons. They are the lightest 
hadrons copiously produced in high energy collisions. 
In the present letter we argue that the pion number 
fluctuations may give a prominent signal of approach- 
ing the BEC point. In fact, there is the BEC line in a 
plane of pion density and temperature. The pion sys- 
tem should be in a metastable state (in thermal, but 
not chemical equilibrium) to reach the BEC line. This 
can be achieved by selecting the samples of events with 
high pion multiplicities. Multipion states are formed 
in high energy nucleus-nucleus collisions, as well as 
in the elementary particle ones. There were several 
suggestions to search for BEC of 7r-mesons (see, e.g., 
Ref. Q). However, complete statistical mechanics cal- 
culations of pion number fluctuations have never been 
presented. There is a qualitative difference in prop- 
erties of the mean multiplicity and of the scaled vari- 
ance of multiplicity fluctuations in different statistical 
ensembles. The results obtained with grand canon- 
ical ensemble (GCE), canonical ensemble (CE), and 
microcanonical ensemble (MCE) for the mean mul- 
tiplicity approach to each other in the large volume 
limit. This reflects the thermodynamic equivalence of 
the statistical ensembles. Recently it has been found 
0) S @ that corresponding results for the scaled vari- 
ance are different in different ensembles, and this dif- 
ference is preserved in the thermodynamic limit. To 



extract the matter properties from analysis of event- 
by-event fluctuations, one needs to fix the samples of 
high energy events, and choose the corresponding sta- 
tistical ensemble for their analysis. This is discussed 
below. 

Let us start with a well known example of non- 
relativistic ideal Bose gas. The occupation numbers, 
n p , of single quantum states, labelled by 3-momenta 
p, are equal to n p = 0,1,... , oo. In the GCE their 
average values, fluctuations, and correlations are the 
following [3]: 

, , 1 



exp 



£r ~ M /T 



- 1 



((A, p f) = W(l + W)^J 
(AripAnk) = v p 5 pk , 



(1) 



where An r 



m denotes the particle mass, 



T and \x are the system temperature and chemical 
potential, respectively (throughout the paper we use 
the units with h — c — k — 1). The average number 
of particles in the GCE reads Q: 
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where V is the system volume. We consider particles 
with spin equal to zero, thus the degeneracy factor 
equals 1. In the thermodynamic limit, V — * oo, the 
sum over momentum states is transformed into the 
momentum integral, J2 P ■ ■ ■ — (V/^ 2 ) . . .p 2 dp. 
This substitution, assumed in all formulae below, is 
valid if the chemical potential in the non-relativistic 
Bose gas is restricted as /i < (or [i < m in relativis- 
tic formulation). When the temperature T decreases 
at fixed ratio, (N)/V, the chemical potential [i in- 
creases and becomes equal to zero at T = Tc, known 
as the BEC temperature. At this point from Eq. $2$ 
one finds, N(V,T=T C , fi=0) = V[mT c /(2TT)} 3 / 2 ((3/2), 
where C(3/2) = 2.612 is the Riemann zeta-function. 
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At p = and T < T c , a macroscopic part, iV (called 
the BE condensate), of the total particle number oc- 
cupies the lowest energy level p = 0. 

Introducing AiV = N — (N) one finds the particle 
number fluctuations in the GCE, 

((AN) 2 ) = ]T (An p An k ) = £ t£ , (3) 
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The limit -/i/T > 1 gives (n p ) < 1. This 
corresponds to the Boltzmann approximation, and 
then from Eqs. (J2EJ) it follows: N(V,T,p) S 
V ex.p(p/T)(mT/2ir) 3 ^ 2 and a; = 1. When // increases 
the scaled variance to becomes larger, u> > 1. This is 
the well known Bose enhancement effect for the par- 
ticle number fluctuations. From Eq. (TJJ at p — ► 
one finds cj — * oo. Two comments are appropri- 
ate. First, for finite systems u> remains finite, and 
u) = oo emerges from Eq. (d]) at // = in the ther- 
modynamic limit — > oo, when the sums over p are 
transformed into the momentum integrals. Second, 
the anomalous fluctuations of the particle number at 
the BEC point correspond to the GCE description. 
In the CE and MCE, the number of particles N in 
a non-relativistic system is fixed by definition, thus, 
ci> c . e , = LUm.c.e. = 0. At p = and T < Tc the average 
number of particles in the BE condensate is propor- 
tional to 1 — (T/Tc) 3 / 2 and it remains the same in the 
CE and MCE at V -> oo. The fluctuations of N are, 
however, very different and strongly suppressed in the 
CE and MCE [|. 

We consider now the relativistic ideal gas of pions, 
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where index j enumerates 3 isospin pion states, 
7r + ,7r - , and 7T°, the energy of one-particle states 
is taken as, e p = (p 2 + to 2 ) 1 / 2 with = 140 
MeV being the pion mass (we neglect a small dif- 
ference between the masses of charged and neutral 
pions). The inequality pj < is a general restric- 
tion in the relativistic Bose gas, and pj = corre- 
sponds to the BEC. In Ref. [5| we discussed in details 
the Bose gas with one conserved charge in the CE 
(V,T,Q — const), i.e. the 7r + 7r~-gas with fixed elec- 
tric charge. This corresponds to the GCE (V,T,/j,q), 
thus, in Eq. §5§ p+ = pq and /i_ = —fiQ for 7r + 
and ir~ , respectively. Approaching the BEC of ir + 
at /iq — > m,r, one finds the relation between Tc and 
Pq = p + — p- (the picture of BEC of ir~ at Q < 
and uq — » —m^ is obtained by a mirror reflection). At 
Tc/nin 1 it coincides with the non-relativistic for- 
mula, Tc = 3.31 Pq m^ 1 , and at Tc/m^ 3> 1 it gives 
T c = 1.73p Q /2 m~ 1/2 (see, e.g., Ref. @). The scaled 
variance uj+ = ((AN + ) 2 ) / (N + ) in the GCE goes to 
infinity. This is similar to the non-relativistic case. 
On the other hand, the scaled variance for negative 



particles, w~ = ((A/V_) 2 )/(/V_), remains finite and 
even decreases with pq. The pion numbers 7V+ and 
7V_ fluctuate in the both GCE and CE. However, the 
exact conservation imposed in the CE on the system 
charge, Q = N + — N-, suppresses anomalous fluctua- 
tions at the BEC point: uj+ e is finite with the upper 
limit, C(2)/C(3) — 1-368, reached at Tc/m — > oo (see 
details in Ref. [f|). 

In what follows we discuss a rather different pion 
system. In the MCE (V, E, Q = 0, iV OT = const) for- 
mulation, the total system energy E, electric charge 
Q = N + — 7V_, and total number of pions N n = 
Nq+N + + N- will be fixed. This corresponds to the 
GCE (V, T, p>Q= 0, p n ) description with p + = p^+PQ, 
P- = p n — pq, and po = p^ in Eq. (JSJ). We restrict 
Pq= and consider BEC when p^-^m^. The pq=0 
corresponds to zero electric charge, Q= or N + = N- , 
in the pion system. 

The pion density is equal to p„(T,p v ) — 
E P j( n p,j)/V- The phase diagram of the ideal pion 
gas in pn — T plane is presented in Fig. [TJ BEC starts 
at T = Tc when p„ — p™ ax — . It gives: 
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where K% is the modified Hankel function. The 
Eq. ([6|) gives the BEC line shown by the solid line 
in Fig. [TJ If T c /m n < 1, from Eq. © one finds, 

T c = 2 7 r[3C(3/2)]- 2 / 3 / 9 2 r /3 TO~ 1 . This corresponds to 
the non-relativistic limit discussed above, but with 
a degeneracy factor g n ~ 3. In the ultrarelativistic 
limit, Tc/mir > 1, from Eq. © it follows: T c = 
[7r 2 /3C(3)] 1 / 3 py 3 . We consider the region in p„ - T 
plane between the p v = and p^ =m w lines. The lines 
of fixed energy density, s(T,p w ) = J2 P ,j e p ( n P,j}/ v > 
are shown as dotted lines in Fig.[T]inside this region for 
three fixed values of e. An increase of p^ at constant 
e leads to the increase of p n and decrease of T. In this 
letter we discuss how the system approaches the BEC 
line — m^, T = Tc), and do not touch the region 
(p v =m 7r ,T<Tc') below this line where the non-zero 
BE condensate is formed. The GCE (V,T, pq, p„), 
MCE (V,E,Q,N n ), and CE (V, T, Q, N„) are equiv- 
alent for average quantities, including average parti- 
cle multiplicities, in the thermodynamic limit. Thus, 
Eq. ([6]) and phase diagram in Fig. 1 remain the same 
in all statistical ensembles. However, the pion number 
fluctuations are very different in different ensembles. 

As an example we consider the high multiplic- 
ity events in p + p collisions at IHEP (Protvino) 
accelerator with the beam energy of 70 GeV (see 
Ref. [10] on the experimental project "Thermaliza- 
tion", team leader V.A. Nikitin). In the reaction 
p + p — * p + p + N n with small final proton mo- 
menta in the c.m.s., the total cm. energy of cre- 
ated pions is E = ^s — 2m p = 9.7 GeV. The trigger 
system designed at JINR (Dubna) selects the events 
with N n > 20 in this reaction. This makes it pos- 
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FIG. 1: The phase diagram of the pion gas with \xq = 0. 
The dashed line corresponds to p 7r (T,[i 7r = 0), and the 
solid line to BEC ([6]). The dotted lines show the states 
with fixed energy densities: e = 6,20,60 MeV/fm 3 . The 
N n numbers in the figure correspond to fi^ = and fi^ — 
at these energy densities for the total pion energy, 
E = 9.7 GeV. 

sible to accumulate the samples of events with fixed 
Nir = 30 -j- 50 and the full pion identification during 
the next 2 years Note that for this reaction the 
kinematic limit is N™ ax = E/m,, 70. The pion 
system in the thermal equilibrium is expected to be 
formed for high multiplicities. The volume of the pion 
gas system is estimated as, V — E/e(T,fj^), and the 
number of pions equals to N„ = Vp w (T, The 
values of N„ at fj^ = and \x v — for 3 different 
values of energy density s are shown in Fig. [TJ 

For Q — 0, the average pion multiplicities, (No) — 
(N±) = Nn/3, are the same in all statistical ensem- 
bles for large systems. This thermodynamic equiva- 
lence is not, however, valid for the scaled variances of 
pion fluctuations. The system with the fixed electric 
charge, Q = 0, the total pion number, N n , and total 
energy of the pion system, E, should be treated in the 
MCE. The volume V is one more (and unknown) MCE 
parameter. The calculations below are carried out in 
a large volume limit, thus, parameter V docs not enter 
explicitly in the formulae for the scaled variances. The 
microscopic correlators for the MCE (V, E, Q — 0, N„) 
equal to: 
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term in the r.h.s. of Eq. corresponds to the GCE. 
Correlations between differently charged pions, j =/=i, 
and between different single modes, p^k, are absent 
in the GCE. It then follows: lo + — lu~ = ui° = to = 
1 + E P (™p) 2 /Ep(«-p) , similar to Eq. (J), but with 
(n p ) = {exp[(Vp 2 +ml - ^)/T] - l}" 1 . In the 




FIG. 2: The scaled variance of neutral pions in the MCE 
is presented as the function of the total number of pions. 
Three solid lines correspond to different energy densities: 
e = 6, 20, 60 MeV/fm 3 . The total energy of the pion sys- 
tem is assumed to be fixed, E — 9.7 GeV. The vertical 
dotted lines correspond to the points on the BEC line at 
the specific values of the energy density. 

GCE the numbers N + , iV_ , and N fluctuate indepen- 
dently of each other. The Bose effects in the pion sys- 
tem are small if p,^ = 0. One finds u) = 1.01 -j- 1.12 in 
the temperature interval T = 40-^160 MeV (note that 
u> = 1 in the Boltzmann approximation). The Bose 
effects increase with fj, w , and lu — > oo at — > m v , i.e. 
approaching the BEC line the GCE calculations give 
anomalous fluctuations for N + , N_, and iVo- The 
MCE (V, E, Q — 0, — const) formulation means 
the restrictions of the fixed total system energy E, 
electric charge Q = N + — 7V_ , and total number of pi- 
ons Nn — Nq + N + + N- for the each microscopic 
state of the system. This changes the pion num- 
ber fluctuations. From Eq. |J7J| one notices that the 
MCE fluctuations of each mode p are reduced, and 
the (anti)correlations between different modes p ^ k 
and between different charge states appear. This re- 
sults in a suppression of all scaled variances u> m .c.e. 
in comparison with the corresponding ones w in the 
GCE. Note that the MCE microscopic correlators Q, 
although being different from that in the GCE, are 
expressed with the quantities calculated in the GCE. 
The MCE scaled variances depend on two GCE pa- 
rameters: T and p^. The straightforward calculations 
lead to the following MCE scaled variance: 
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Due to conditions, iV + = iV_ and N + +N-+N Q = N v 



it follows, w± c e = w^, 
where N ch = N + + N- 



/4 and w. 
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The behavior of w^. c . e . is 



shown in Fig. [5] To make a correspondence with N n 
values, we consider again the p + p — > p + p + N n colli- 
sions at the beam energy of 70 GeV and take the pion 
system energy to be equal to E = 9.7 GeV. Despite 
of the MCE suppression the scaled variances for the 
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number fluctuations of ir° and ir increase dramati- 
cally and abruptly when the system approaches the 
BECline. 

As an instructive example let us consider the MCE 
(V, E, Q— 0, N c h — const), i.e. fixed N c h instead of 
N^. The corresponding GCE formulation gives the 
following pion chemical potential: u + = fj, n , /xq = 0, 
//_ = /ijr in Eq. © ((1q = 0, as before, because of 
Q = condition). When — > the system ap- 
proaches the BEC line for ir + and ir~ . The thermo- 
dynamic behavior and position of this BEC line can be 
easily found. Approaching the BEC line one can also 
find u>~^ ce — ► oo. The pion number fluctuations are, 
however, very different in the MCE (V, E, Q=0, N ch ). 
In the statistical ensembles with fixed N c h and Q no 
anomalous BEC fluctuations are possible. The num- 
bers of N + and 7V_ are completely fixed by the condi- 
tions Q = N + -N- =0 and N ch = N+ + N- = const, 
thus, = ui^ c e — 0. The number N fluctuates, 
but fiQ = 0, thus, neutral pions are far away from the 
BEC line and their fluctuations are small, oj° w 1 in 
all statistical ensemble formulations. 

The broad distributions over N and N c h close to 
the BEC line also implies large fluctuations of the 
/ = No/ N c h ratio. These large fluctuations were sug- 
gested (see, e.g., Ref. [Hj]) as a possible signal for the 
disoriented chiral condensate (DCC). The DCC leads 
to the distribution of / in the form, dW(f)/df = 
l/(2v / 7). The thermal Bose gas corresponds to the 
/-distribution centered at / = 1/3. Therefore, /- 
distributions from BEC and DCC are very different, 
and this gives a possibility to distinguish between 
these two phenomena. 

The calculations presented in this letter should be 
improved by taking into account the finite size ef- 
fects, pion-pion interactions, and some other effects. 
However, the described BEC scenario may survive the 
complications. A crucial point is the analysis of the 



samples of high N„ events. The following inequalities 
are always hold for particle number fluctuations in dif- 
ferent ensembles: w % mce < w* >e< < uj % c er Therefore, 
if the anomalous BEC fluctuations are present in the 
MCE, they are also exist (and even larger) in the CE 
and GCE. The reverse statement is not true. The 
anomalous BEC fluctuations of the GCE may disap- 
pear in the CE or MCE. We found that for the system 
with N^ = const and Q = the anomalous BEC fluc- 
tuations do not wash out by exact conservation laws 
of the CE and MCE. The required N n values for the 
BEC are much larger than the average pion multiplic- 
ity per collision, thus, these high N n events are rather 
rare and give negligible contributions to inclusive ob- 
servables in high energy collisions. With increasing 
of Nn in the sample with fixed total energy, the tem- 
perature of the pion system has to decrease and it 
approaches the BEC line. This can happen in differ- 
ent ways: at constant energy density e, at constant 
pion density p v , or with decreasing of both e and p n . 
The pion system should move to the BEC line one 
way or another. In the vicinity of the BEC line (no 
BE condensate is yet formed) one observes an abrupt 
and anomalous increase of the scaled variances of neu- 
tral and charged pion number fluctuations. This could 
(may be even should) be checked experimentally. 
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